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Effect of a static phase transition on searching dynamics 
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We consider a one dimensional Euclidean network which is grown using a preferential attachment. 
Here the jth incoming node gets attached to the ith existing node with the probability Hi oc kilfj, 
where Uj is the Euclidean distance between them and ki the degree of the ith node. This network is 
known to have a static phase transition point at a c ~ 0.5. On this network, we employ three different 
searching strategies based on degrees or distances or both, where the possibility of termination of 
search chains is allowed. A detailed analysis shows that these strategies are significantly affected by 
the presence of the static critical point. The distributions of the search path lengths and the success 
rates are also estimated and compared for the different strategies. These distributions appear to be 
marginally affected by the static phase transition. 
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I. INTRODUCTION 

Static critical points are known to affect dynamical 
phemomena usually giving rise to critical slowing down, 
e.g., as in the relaxation of the order parameter in mag- 
netic systems. Dynamical critical phenomena is a well 
studied and important topic, exploring the dynamical 
behaviour of systems, especially at the thermal critical 
point p|. In many systems, phase transitions driven by 
factors other than temperature can occur, as for example 
the geometrical phase transition occurring in percolation. 
Even in magnetic systems, e.g., in the axial next near- 
est neighbour Ising (ANNNI) model, a phase transition 
occurs at zero termperature when the competing second 
neighbour interaction takes up a certain value. How- 
ever, here a zero temperature quenching dynamics fails 
to carry any signature of the phase transition 0] . 

Recently, with the discovery of the small world effect 
in real networks, many theoretical models have been set 
up to mimic small worlds. In some of these models, in- 
teresting phase transitions have been noted, e.g., in the 
Watts-Strogatz model [1 E| > where one starts with near- 
est neighbour links only and then rewires links to long 
range neighbours with probability p, the small world ef- 
fect is observed even as p — > 0. Phase transitions in 
models in which the linking probability is dependent on 
spatial and/or temporal factors have also been observed 
d There is no temperature associated with these 
networks. Of course, if one considers spin systems, e.g., 
the Ising model, on such networks, it is possible to ob- 
tain thermal phase transitions as well. Dynamical studies 
of such systems at both zero and non-zero temperature 
have shown unexpected phenomena, as for example freez- 
ing in case of the quenching dynamics of the Ising model 
on small world networks [7j. 

The idea of a small world first emerged from a real 
dynamical experiment made on the US population by 
Milgram [8| which showed that on an average there are 
six steps required to reach another individual. Later, 



a mathematical definition of small world property was 
proposed; by the small world property it is meant that if 
any two nodes in the network is separated by an average 
number of s steps, then s oc ln(AT), where N is the total 
number of nodes in the network. In some networks, even 
slower variation (i.e., sub-logarithmic) scaling has been 
observed @. 

Following Milgram's original experiment, several new 
experiments have been done to study the searching dy- 
namics in real social networks pi Qj|- A considerable 
number of theoretical studies on searching phenomena 
has also been made recently pj, El EI El 111 E3, El 

While the small world property is a "static" property, 
i.e., calculated on the basis of global knowledge and with- 
out any dynamics involved, it must be noted that it is not 
necessary that a navigation or search on a small world 
network would show the small world property, i.e., the 
dynamic paths Sd may not scale as ln(AT). This is be- 
cause searching is done using local information only. This 
was explicitly shown by Kleinberg [15[ in a theoretical 
study where nodes where placed on a two dimensional 
Euclidean space. Each node here has connections to its 
nearest neighbours as well as to neighbours at a distance 
I with probability 



P(l) oc r 



(1) 



Although the network is globally a small world for a 
range of values of a, navigation on such networks using 
greedy algorithm showed a small world behaviour with 
Sd ~ (log(N)) 2 only at a = 2. In general the path length 
Sd showed a sublinear power law increase with N. 

Although searching dynamics is not comparable to dy- 
namics of relaxation, quenching etc., still we have some 
indication that it bears the signature of a static phase 
transition from some earlier studies pi E3 • ^ n one di- 
mension, the Euclidean network in which connections are 
made using ([T]) has been shown to have three phases cor- 
responding to the scaling behaviour of the link lengths. It 
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is a small world for for a < 1 where the link lengths scale 
as N; finite dimensional network for 1 < a < 2 where the 
link lengths scale sublinearly with N and a regular net- 
work for a > 2 where the link lengths are O(l) [2a]. Here 
the existence of the two phase transitions are quite con- 
firmed and the greedy searching strategy has been shown 
to bear the signatures of both [18| ■ Correspondingly, it is 
expected that in two dimensions, there is a transition to 
the small world behaviour at a c = 2 and the result that 
the search path lengths scale uniquely at this value of 
a shows that the searching dynamics is indeed sensitive 
to it. However, this has been observed for a particular 
algorithm and may not be true always. In this paper, 
we have considered a network where different kind of al- 
gorithms can be used to check the sensitivity to a static 
phase transition. 

We have considered realistic searches on this network. 
As in 24], here also the search paths may terminate. 
Thus there is a success rate also involved in the study. 
To study the effect of the parameter a which governs 
the phase transition in the model (details given in the 
next section), we have computed the path length, the 
success rate and the ratio of the two as well. The last 
quantity shows a power law variation with the network 
size, N~ s , and can be taken as a reliable measure to 
compare different search strategies [HI ■ We find out that 
the variation of 5 with a indeed shows the signature of the 
phase transition to a considerable extent for the different 
algorithms. 

We describe the model and the algorithms in section II 
and the results in section III. We have computed the dis- 
tribution of the path lengths and success rates for all the 
algorithms also which are presented in sec IV. Summary 
and concluding remarks are made in sec V. 

II. THE NETWORK AND THE STRATEGIES 

We have considered a growing Euclidean network in 
which the nodes are added one by one using a preferen- 
tial attachment such that the probability that an incom- 
ing node j gets attached to an existing node i with ki 
neighbours at that time is [f| 

Hi cx kl~ a , (2) 

where Uj is the distance between the ith and the jth 
nodes. The nodes are generated randomly at sites 
Xi (0 < xi < 1) on a continuous one dimensional line and 
each new node gets attached to m pre-existing nodes. 

For m = 1, this network was shown to have a phase 
transition at a c : below a c it has a stretched exponential 
degree distribution while above it the degree distribution 
has a power law tail with the exponent 7 equal to 3.0. 
The value of a c was found to be close to 0.5 [fl- 
it is to be noted that this model, in the limit a = 
is nothing but the scale-free Barabasi Albert model [26| . 
On scale-free networks, a degree based search is natu- 
ral to adopt while the greedy algorithm appears to be 



the most popular one from the findin gs o f the original 
Milgram experiment as well as that of [id fllj on a Eu- 
clidean network. We have therefore employed three dif- 
ferent search algorithms here, based on the degrees of the 
nodes, or distances, or both. 

All the search strategies follow the basic rules : 
After a source node and a target node are selected ran- 
domly, the source node will send the signal to one of its 
neighbouring nodes provided that node has not already 
taken part in the search. This is repeated till the message 
reaches a node which is connected to the target node and 
this scheme is in tune with Milgram-like experiments. In 
course of this search, it may happen that a node cannot 
pass the signal to any of its neighbour as they have al- 
ready taken part in the search. In that case, the search 
terminates at that node. Such searches are termed un- 
successful. The fraction of successful searches is called 
the success rate p. The average number of steps taken 
to reach the target in a successful search is the average 
dynamic path length Sd- We also calculate the quantity 
u, defined as « = -. 

The choice of the neighbour to whom the signal is 
being passed depends on the strategy. The three search 
strategies considered in the present work are as follows: 

(1) Highest Degree Search (HDS): Here after a source 
and a target pair are chosen randomly the source scans its 
nearest neighbours and chooses the one with the highest 
degree to pass on the signal. 

(2) Nearest neighbour search (NNS): In this strat- 
egy, after the source-target pair is chosen randomly, the 
source chooses from among its nearest neighbours, the 
one whose Euclidean distance (I) from the target is the 
least. It may be noted that in conventional greedy algo- 
rithms, the strategy is to pass the message to a neigh- 
bour which is nearer to the target than itself. In the 
present case, this condition has not been imposed and 
therefore in an intermediate stage, the message may "pro- 
ceed backward" . This is analogous to allowing a system 
to go a higher energy configuration in simulated anneal- 
ing applied to the dynamics of frustrated systems like 
spin glasses and to combinatorial optimisation problems. 

(3) Optimised Search (OS): In this strategy, we follow 
an algorithm where the degree (fc) of a node as well as 
its Euclidean distance (/) from the target are taken into 
account. Here, after a source-target pair is picked up at 
random, the ratio £ = k/l is calculated for all the nearest 
neighbours of the source and the one with the highest 
value of £ is chosen to pass on the signal. 



III. RESULTS FOR p, s d AND 8 

We have simulated the networks with a maximum of 
N = 5000 nodes using upto 1000 configurations. For each 
configuration, the searching is repeated N/2 times with 
randomly chosen source-target pairs. We have considered 
two cases, m = 1 and m = 2. 
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FIG. 1: Variation of p, Sd and p with a for the three search 
strategies for N = 1000. The parameter a is varied from 
— 10.0 to +10.0. A transitional behaviour is observed for all 
the strategies around the static critical point of the system, 
i.e., near a = a c ~ 0.5 



A. Case I, m = 1 

Here the tunable parameter a has been varied from 
— 10.0 to +10.0. Once the network is generated following 
eq. @, and the navigation has been simulated following 
one of the three strategies described in the last section, 
the success rate p and the average search length Sd are 
evaluated and their variation with a and N is noted. 



First we have made a comparison of the three strate- 
gies by analysing the variation of p, Sd and \i with a 
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FIG. 2: Variations of p with system size N are shown for 
a — 1.0 for the three search strategies for m = 1 and m — 2. 
The success rate decreases with increasing system size. 
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FIG. 3: Variations of the average path length Sd with system 
size iV for the three different strategies are shown for a — 1.0 
for m = 1 and m = 2. It is observed that the average search 
length increases with increasing N. 



for a fixed system size N. 
variations for N — 1000. 



We show in Fig. 1 these 



For the HDS strategy, it is observed that p shows a 
slow increase as a increases from —10.0 upto 0.0 after 
which it falls sharply till a = 2.0 and finally tends to 
saturate beyond a = 5.0. The value of Sd for this strategy 
however remains constant from a = —10.0 to a = 0.0 and 
increases slowly from this value also showing a tendency 
to saturate at large values of a. 

For the NNS strategy, p and Sd remain very small for 
a < 0; p shows a gradual increase between a — 0.0 and 
a = 2.0. The values of Sd however increase quite rapidly 
between a = 0.0 and a = 2.0. 

For the OS strategy, p has a slow decrease from a = 
— 10.0 upto a — 0.0 and then it increases quite sharply 
between a = 0.0 and a = 5.0 beyond which it saturates. 
Similarly Sd decreases slowly between a = —10.0 and 
a = 0.0, then increases very sharply upto a = 2.0. 

A saturation of both p and Sd is expected for all the 
algorithms as the network approaches the behaviour 
of a growing network in which links are made to the 
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FIG. 4: The ratio p = p/sd vs N plots at a = 1.0 for m = 1 
and m = 2. 



nearest neighbours (i.e., the a — > oo limit) for large a. 
Similar saturation behaviour for the static properties 
was observed in [j| 

Looking at Sd alone, it would seem that the HDS is still 
the best strategy even at a » 1, when the network is 
not scale-free. However, p for HDS becomes very low here 
indicating that very few chains are completed, in which 
case chains tend to be 'short'. This explains the above 
observation for Sd- Similarly, for a < 0, NNS would seem 
to be the best from the values of s<j. On the other hand, 
from the p plots, OS seems best for a >> 1, while HDS 
seems best for a < 0. As in (24|, here also we compute 
p, which incorporates both p and Sd, to comment on the 
relative capabilities of the three strategies. From the p, 
vs a plots above it is apparent that HDS works best upto 
a ~ 0.5 while OS is best for a > 0.5 for this particular 
value of N. NNS works rather poorly for a < and 
performs relatively better for a > 0. 

We find that the behaviour of p in general closely fol- 
lows that of the success rates p. This may indicate that 
rather than the path lengths, which are "small" in all 
cases, the success rate decides the quality of the search 
strategy here. 

Next we discuss the behaviour of the above quantities 
with N. 

We show typical plots of p and Sd against N (Figs 2,3) 
for a fixed value of a = 1.0. p clearly shows a power law 
decay with N . Sd apparently has a power law increase, 
with a very small exponent (~ 0.01). However here one 
expects Sd ~ ln(iV) as the network has a tree structure 
for m = 1. Indeed, we find that the exponent tends to 
decrease at larger N, consistent with this. 




FIG. 5: 

Comparison of the different search strategies showing the 
variation of the exponent S with a. All the three search 
strategies show transitional behaviour close to the static 
phase transition point a ~ 0.5. m = 1 here. 

From Fig. 4, we find that p, shows a power law varia- 
tion with the system size N, 

-6 



(3) 



where the value of S varies with a. 

We have computed 6 for different values of a (from 
a = —10.0 to a = +10.0) and shown its variation against 
a for all the three strategies in Fig. 5. A smaller value 
of S indicates a more successful strategy. For values of 
a < 0.5, i.e., when the network is still scale free, the 
strategy purely dependent on degree works better com- 
pared to those dependent on distance. For a >> on 
the other hand, the distance dependent searches perform 
better than the purely degree dependent search. In fact, 
the OS appears to be the best strategy immediately be- 
yond the static critical point while the NNS works better 
than the HDS only when a > 1.5. For a < a c the net- 
work is scale free and there are several high degree nodes 
present so that the HDS strategy wins over the other 
two. However this strategy becomes inefficient beyond 
a > a c , when the system is no longer scale free and high 
degree nodes are no longer available. On the other hand, 
in this region, distance-based search strategies work more 
efficiently as the network has nodes linked to closer neigh- 
bours and both the algorithms, NNS and OS are greedy 
algorithms as far as distances are concerned. Although 
for a > a c , the OS strategy works best, the exponent 5 is 
never very close to zero, which means that the dynamic 
small world effect [24| is absent here. For all the strate- 
gies however, < S < 1, consistent with the boundary 
values obtained in [24]. It is observed that for a nar- 
row region of values of a > a c , HDS is still better than 
NNS which indicates that the relevance of the degree of a 
node reduces gradually, once the network becomes non- 
scale free. The fact that the OS performs best even for 
very large values of a also suggests that the degree is 
never totally irrelevant. 

For both HDS and OS, S shows a drastic in- 
crease/decrease, indicating a sharp transition at a ~ 0.5, 
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which is the static phase transition point. Unlike HDS 
and OS, there is no sharp change in behaviour in 5 for 
NNS and it is affected by the static phase transition point 
to a lesser extent compared to the other two strategies. 

NNS at a — is nothing but a random search. 6 
for NNS remains almost a constant for — oo < a < 0, 
showing that it is never better than a random search 
strategy for this region. For negative values of a, nodes at 
large distance are linked up, but it does not help a greedy 
algorithm. Incorporating k in the algorithm surely helps 
as OS is better than NNS here. In fact S shows a variation 
with a for a < only for the OS strategy. For HDS, the 
plot of 6 versus a is close to a perfect sigmoid, showing 
accountable variation only around a = a c . 



B. Case II, m = 2 

As long as m = 1, the network cannot have any loop 
and the path from one node to another is unique. To 
introduce loops to the lowest order, we have next consid- 
ered searches on networks generated using ^ once again 
where each incoming node can get two links (m = 2). It 
is expected that the static phase transition point remains 
same for m = 2. 

With loops, the success rate should be higher but the 
search lengths may increase several times. In the last 
subsection, we found that the HDS and OS are the more 
effective strategies and we have used only these two in 
the present study. The variation of p, Sd and p with 
N for m = 2 have been plotted in figures 2, 3 and 4 
along with the m — 1 plots. All these show power law 
variations with N . As expected, we find a slower decay of 
p compared to that for m = 1 while Sd increases clearly 
with a power law compared to the logarithmic increase 
obtained for m = 1. The results show that 5 for a = 0.0 
and a = — 1 are very close for OS and HDS; in fact for 
both these values of a, S ~ 0.7 for the two strategies. 
This is greater than that of the m = 1 case. A higher 
value of 5 indicates a deterioration in performance, the 
rapid increase in the path lengths being the reason behind 
this deterioration. 

For a — 1 on the other hand, when the network is no 
longer a scale-free network, 5 values are drastically dif- 
ferent for the two strategies. For OS, it is much smaller, 
~ 0.3, while for HDS, it is around 0.8. The reason be- 
hind this is, with the OS strategy, the path lengths scale 
with a much smaller exponent in the non scale- free region 
(Fig. 3). 



IV. DISTRIBUTIONS 

In the last section, we have compared some search 
strategies on a growing Euclidean network by comput- 
ing the quantities like success rates and search lengths, 
where the mean value of these quantities have been used 
to obtain the scaling behaviour. We have also computed 
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FIG. 6: Plots of the distributions of p for the three strategies 
for a = 0.0, 5.0 and —5.0 are shown for m = 1. The system 
size used is N = 2000. 



the distributions of these quantities to see whether the 
presence of the phase transition has any effect on these. 
The results for three values of a, a < 0, a — and 
a > a c are reported here. 

The distributions of the success rate p for the three 
strategies with m = 1 show the following general fea- 
tures (Fig. 6) 

1. All of them have a well defined peak. 

2. They are symmetric. 

3. There is no long tail. 

4. Distributions are skewed when the mean value is small, 
in fact very few points with non-zero value appear here. 
However, when the mean value is larger, there is a suffi- 
cient broadening, no matter which strategy is being used. 

Overall, we do not find any indication that the static 
phase transition point has a significant influence on the 
form of the distributions. Since the data points are few, 
we do not attempt a fitting but in all probability these 
distributions are gaussian or nearly gaussian. 

The distribution of the path length Sd, on the other- 
hand, is definitely not symmetric (Fig. 7) for any of the 
strategies at any value of a. For m = 1, it has a broad 
peak. None of the distributions have a long tail. For 
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FIG. 7: The distributions of Sd for the three strategies for 
a = 0.0, 2.0 and —2.0 are shown in a log-linear scale. The 
system size used is N = 2000. m = 1 here. 




50 100 150 200 250 300 350 




50 100 150 200 250 300 350 400 




50 100 150 200 250 300 



FIG. 8: The distributions of average search lengths for OS 
and HDS strategy at a — 1,0,-1 when m = 2 are shown. 
N = 2000 here. 



m = 2, when the success rate becomes much higher, we 
find the peaks shifted towards lesser values of Sd (Fig. 8), 
consistent with the observation of [3]. The presence of 
a peak at smaller values of Sd for both m = 1 and m = 2 
also show that shorter paths are more probable [ll| . 

Again, for m = 1, the number of data points are few 
and larger fluctuations exist making it difficult to fit the 
data to any familiar functional form. For m = 2, there 
is a larger number of points and one can immediately 
see that the optimised search strategy has a clear-cut 
exponential decay when a > a c whereas for a < a c , it 
has a slower than exponential decay. However, no such 
change in behaviour is observed for the HDS strategy, it 
is slower than exponential in each case. 



V. SUMMARY AND CONCLUSIONS 

In this work, we have applied different search strate- 
gies to a network which undergoes a phase transition from 
a scale-free to a non scale-free phase. One of our aims 
was to investigate whether such a phase transition signifi- 
cantly affects the search or not as in purely Euclidean net- 



works, such an indication is there [T^, [ll|. We find that 
the search strategies indeed show a significant change in 
behaviour at or near the phase transition to different ex- 
tents. It appears that the degree based searches are more 
sensitive to the phase transition. 

The searching scheme used here allows termination 
of messages as the restriction that the message can be 
passed only once by any messenger has been imposed. In 
reality, of course, several other reasons may exist for a 
termination [ll[ . The analysis of the results therfore has 
been made based on an approach recently suggested by 
one of us 24], in which both search paths and success 
rates are taken into consideration. 

Searching phenomena is vastly studied in social net- 
works and the present study also uses ideas relevant to 
Milgram-like searches. Most social networks being non 
scale-free, our results for the network in its non scale- 
free phase is important in the context of social searches. 
Here the best performance is shown by the optimised 
search (OS) strategy, in which a node sends the signal 
to its neighbour having the largest value of k/l (k is its 
degree and I the distance from the target, see sec II). We 
have used here three strategies and for none of them we 
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observe a dynamical small world effect, i.e., S is never 
very close to zero. However, our list of strategies is by no 
means exhaustive. The OS scheme can be generalised by 
making the message passing rule that a node sends the 
signal to a neighbour with the largest value of k a /l b , in- 
troducing tunable parameters a and b. In this paper, we 
have only considered the limiting cases a —> oo (HDS), 
b -> oo (NNS) and a = 1, b = 1 (OS). It may be an in- 
teresting future study to find out whether in the a — b 
plane, one obtains regions of dynamic small world effect. 
A more detailed study for m — 2 (or more) can also be 



done for which we have presented results at some specific 
values of a only. 

We have also estimated the distributions for the suc- 
cess rate and path lengths. The static phase transition 
seems to seriously affect only the distribution for the path 
lengths for the optimised strategy when m = 2. 
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